DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues.

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Expanding the PRAAS spectrum: De novo mutations of immunoproteasome subunit ß-type 10 in six infants with SCID-Omenn syndrome.

Mutations in proteasome ß-subunits or their chaperone and regulatory proteins are associated with proteasome-associated autoinflammatory disorders (PRAAS). We studied six unrelated infants with three de novo heterozygous missense variants in PSMB10, encoding the proteasome ß2i-subunit. Individuals presented with T-B-NK± severe combined immunodeficiency (SCID) and clinical features suggestive of Omenn syndrome, including diarrhea, alopecia, and desquamating erythematous rash. Remaining T cells had limited T cell receptor repertoires, a skewed memory phenotype, and an elevated CD4/CD8 ratio. Bone marrow examination indicated severely impaired B cell maturation with limited V(D)J recombination. All infants received an allogeneic stem cell transplant and exhibited a variety of severe inflammatory complications thereafter, with 2 peri-transplant and 2 delayed deaths. The single long-term transplant survivor showed evidence for genetic rescue through revertant mosaicism overlapping the affected PSMB10 locus. The identified variants (c.166G>C [p.Asp56His] and c.601G>A/c.601G>C [p.Gly201Arg]) were predicted in silico to profoundly disrupt 20S immunoproteasome structure through impaired ß-ring/ß-ring interaction. Our identification of PSMB10 mutations as a cause of SCID-Omenn syndrome reinforces the connection between PRAAS-related diseases and SCID.

PANDORA v2.0: Benchmarking peptide-MHC II models and software improvements.

T-cell specificity to differentiate between self and non-self relies on T-cell receptor (TCR) recognition of peptides presented by the Major Histocompatibility Complex (MHC). Investigations into the three-dimensional (3D) structures of peptide:MHC (pMHC) complexes have provided valuable insights of MHC functions. Given the limited availability of experimental pMHC structures and considerable diversity of peptides and MHC alleles, it calls for the development of efficient and reliable computational approaches for modeling pMHC structures. Here we present an update of PANDORA and the systematic evaluation of its performance in modelling 3D structures of pMHC class II complexes (pMHC-II), which play a key role in the cancer immune response. PANDORA is a modelling software that can build low-energy models in a few minutes by restraining peptide residues inside the MHC-II binding groove. We benchmarked PANDORA on 136 experimentally determined pMHC-II structures covering 44 unique αß chain pairs. Our pipeline achieves a median backbone Ligand-Root Mean Squared Deviation (L-RMSD) of 0.42 Å on the binding core and 0.88 Å on the whole peptide for the benchmark dataset. We incorporated software improvements to make PANDORA a pan-allele framework and improved the user interface and software quality. Its computational efficiency allows enriching the wealth of pMHC binding affinity and mass spectrometry data with 3D models. These models can be used as a starting point for molecular dynamics simulations or structure-boosted deep learning algorithms to identify MHC-binding peptides. PANDORA is available as a Python package through Conda or as a source installation at https://github.com/X-lab-3D/PANDORA.

Understanding structure-guided variant effect predictions using 3D convolutional neural networks.

Predicting pathogenicity of missense variants in molecular diagnostics remains a challenge despite the available wealth of data, such as evolutionary information, and the wealth of tools to integrate that data. We describe DeepRank-Mut, a configurable framework designed to extract and learn from physicochemically relevant features of amino acids surrounding missense variants in 3D space. For each variant, various atomic and residue-level features are extracted from its structural environment, including sequence conservation scores of the surrounding amino acids, and stored in multi-channel 3D voxel grids which are then used to train a 3D convolutional neural network (3D-CNN). The resultant model gives a probabilistic estimate of whether a given input variant is disease-causing or benign. We find that the performance of our 3D-CNN model, on independent test datasets, is comparable to other widely used resources which also combine sequence and structural features. Based on the 10-fold cross-validation experiments, we achieve an average accuracy of 0.77 on the independent test datasets. We discuss the contribution of the variant neighborhood in the model's predictive power, in addition to the impact of individual features on the model's performance. Two key features: evolutionary information of residues in the variant neighborhood and their solvent accessibilities were observed to influence the predictions. We also highlight how predictions are impacted by the underlying disease mechanisms of missense mutations and offer insights into understanding these to improve pathogenicity predictions. Our study presents aspects to take into consideration when adopting deep learning approaches for protein structure-guided pathogenicity predictions.

Parallel evolution and enhanced virulence upon in vivo passage of an RNA virus in Drosophila melanogaster.

Virus evolution is strongly affected by antagonistic co-evolution of virus and host. Host immunity positively selects for viruses that evade the immune response, which in turn may drive counter-adaptations in host immune genes. We investigated how host immune pressure shapes virus populations, using the fruit fly Drosophila melanogaster and its natural pathogen Drosophila C virus (DCV), as a model. We performed an experimental evolution study in which DCV was serially passaged for ten generations in three fly genotypes differing in their antiviral RNAi response: wild-type flies and flies in which the endonuclease gene Dicer-2 was either overexpressed or inactivated. All evolved virus populations replicated more efficiently in vivo and were more virulent than the parental stock. The number of polymorphisms increased in all three host genotypes with passage number, which was most pronounced in Dicer-2 knockout flies. Mutational analysis showed strong parallel evolution, as mutations accumulated in a specific region of the VP3 capsid protein in every lineage in a host genotype-independent manner. The parental tyrosine at position ninety-five of VP3 was substituted with either one of five different amino acids in fourteen out of fifteen lineages. However, no consistent amino acid changes were observed in the viral RNAi suppressor gene 1A, nor elsewhere in the genome in any of the host backgrounds. Our study indicates that the RNAi response restricts the sequence space that can be explored by viral populations. Moreover, our study illustrates how evolution towards higher virulence can be a highly reproducible, yet unpredictable process.

Repurposing Drugs Based on Evolutionary Relationships Between Targets of Approved Drugs and Proteins of Interest.

Drug repurposing has garnered much interest as an effective method for drug development among biopharmaceutical companies. The availability of information on complete sequences of genomes and their associated biological data, genotype-phenotype-disease relationships, and properties of small molecules offers opportunities to explore the repurpose-able potential of existing pharmacopoeia. This method gains further importance, especially, in the context of development of drugs against infectious diseases, some of which pose serious complications due to emergence of drug-resistant pathogens. In this article, we describe computational means to achieve potential repurpose-able drug candidates that may be used against infectious diseases by exploring evolutionary relationships between established targets of FDA-approved drugs and proteins of pathogen of interest.

Exploring anti-malarial potential of FDA approved drugs: an in silico approach.

BACKGROUND: The critically important issue on emergence of drug-resistant malarial parasites is compounded by cross resistance, where resistance to one drug confers resistance to other chemically similar drugs or those that share mode of action. This aspect requires discovery of new anti-malarial compounds or formulation of new combination therapy. The current study attempts to contribute towards accelerating anti-malarial drug development efforts, by exploring the potential of existing FDA-approved drugs to target proteins of Plasmodium falciparum. METHODS: Using comparative sequence and structure analyses, FDA-approved drugs, originally developed against other pathogens, were identified as potential repurpose-able candidates against P. falciparum. The rationale behind the undertaken approach is the likeliness of small molecules to bind to homologous targets. Such a study of evolutionary relationships between established targets and P. falciparum proteins aided in identification of approved drug candidates that can be explored for their anti-malarial potential. RESULTS: Seventy-one FDA-approved drugs were identified that could be repurposed against P. falciparum. A total of 89 potential targets were recognized, of which about 70 are known to participate in parasite housekeeping machinery, protein biosynthesis, metabolic pathways and cell growth and differentiation, which can be prioritized for chemotherapeutic interventions. An additional aspect of prioritization of predicted repurpose-able drugs has been explored on the basis of ability of the drugs to permeate cell membranes, i.e., lipophilicity, since the parasite resides within a parasitophorous vacuole, within the erythrocyte, during the blood stages of infection. Based on this consideration, 46 of 71 FDA-approved drugs have been identified as feasible repurpose-able candidates against P. falciparum, and form a first-line for laboratory investigations. At least five of the drugs identified in the current analysis correspond to existing antibacterial agents already under use as repurposed anti-malarial agents. CONCLUSIONS: The drug-target associations predicted, primarily by taking advantage of evolutionary information, provide a valuable resource of attractive and feasible candidate drugs that can be readily taken through further stages of anti-malarial drug development pipeline.

Comparison of Leptospira interrogans and Leptospira biflexa genomes: analysis of potential leptospiral-host interactions.

Leptospirosis, a potentially life-threatening disease, remains the most widespread zoonosis caused by pathogenic species of Leptospira. The pathogenic spirochaete, Leptospira interrogans, is characterized by its ability to permeate human host tissues rapidly and colonize multiple organs in the host. In spite of the efforts taken to comprehend the pathophysiology of the pathogen and the heterogeneity posed by L. interrogans, the current knowledge on the mechanism of pathogenesis is modest. In an attempt to contribute towards the same, we demonstrate the use of an established structure-based protocol coupled with information on subcellular localization of proteins and their tissue-specificity, in recognizing a set of 49 biologically feasible interactions potentially mediated by proteins of L. interrogans in humans. We have also presented means to adjudge the physicochemical viability of the predicted host-pathogen interactions, for selected cases, in terms of interaction energies and geometric shape complementarity of the interacting proteins. Comparative analyses of proteins of L. interrogans and the saprophytic spirochaete, Leptospira biflexa, and their predicted involvement in interactions with human hosts, aided in underpinning the functional relevance of leptospiral-host protein-protein interactions specific to L. interrogans as well as those specific to L. biflexa. Our study presents characteristics of the pathogenic L. interrogans that are predicted to facilitate its ability to persist in human hosts.

Computational recognition and analysis of hitherto uncharacterized nucleotide cyclase-like proteins in bacteria.

UNLABELLED: Evolutionary relationship between class III nucleotide cyclases and an uncharacterized set of bacterial proteins from Actinobacteria, Bacteroidetes and Proteobacteria has been recognized and analyzed. Detailed analyses of sequence and structural features resulted in the recognition of potential cyclase function conferring residues and presence of signature topological motif (ßααßßαß) in the uncharacterized set of bacterial proteins. Lack of transmembrane domains and signal peptide cleavage sites is suggestive of their cytosolic subcellular localization. Furthermore, analysis on evolutionarily conserved gene clusters of the predicted nucleotide cyclase-like proteins and their evolutionary relationship with nucleotide cyclases suggest their participation in cellular signalling events. Our analyses suggest expansion of class III nucleotide cyclases. REVIEWERS: This article was reviewed by Eugene Koonin and Michael Gromiha.

Recognizing drug targets using evolutionary information: implications for repurposing FDA-approved drugs against Mycobacterium tuberculosis H37Rv.

Drug repurposing to explore target space has been gaining pace over the past decade with the upsurge in the use of systematic approaches for computational drug discovery. Such a cost and time-saving approach gains immense importance for pathogens of special interest, such as Mycobacterium tuberculosis H37Rv. We report a comprehensive approach to repurpose drugs, based on the exploration of evolutionary relationships inferred from the comparative sequence and structural analyses between targets of FDA-approved drugs and the proteins of M. tuberculosis. This approach has facilitated the identification of several polypharmacological drugs that could potentially target unexploited M. tuberculosis proteins. A total of 130 FDA-approved drugs, originally intended against other diseases, could be repurposed against 78 potential targets in M. tuberculosis. Additionally, we have also made an attempt to augment the chemical space by recognizing compounds structurally similar to FDA-approved drugs. For three of the attractive cases we have investigated the probable binding modes of the drugs in their corresponding M. tuberculosis targets by means of structural modelling. Such prospective targets and small molecules could be prioritized for experimental endeavours, and could significantly influence drug-discovery and drug-development programmes for tuberculosis.

SInCRe-structural interactome computational resource for Mycobacterium tuberculosis.

We have developed an integrated database for Mycobacterium tuberculosis H37Rv (Mtb) that collates information on protein sequences, domain assignments, functional annotation and 3D structural information along with protein-protein and protein-small molecule interactions. SInCRe (Structural Interactome Computational Resource) is developed out of CamBan (Cambridge and Bangalore) collaboration. The motivation for development of this database is to provide an integrated platform to allow easily access and interpretation of data and results obtained by all the groups in CamBan in the field of Mtb informatics. In-house algorithms and databases developed independently by various academic groups in CamBan are used to generate Mtb-specific datasets and are integrated in this database to provide a structural dimension to studies on tuberculosis. The SInCRe database readily provides information on identification of functional domains, genome-scale modelling of structures of Mtb proteins and characterization of the small-molecule binding sites within Mtb. The resource also provides structure-based function annotation, information on small-molecule binders including FDA (Food and Drug Administration)-approved drugs, protein-protein interactions (PPIs) and natural compounds that bind to pathogen proteins potentially and result in weakening or elimination of host-pathogen protein-protein interactions. Together they provide prerequisites for identification of off-target binding.

Homology-Based Prediction of Potential Protein-Protein Interactions between Human Erythrocytes and Plasmodium falciparum.

Plasmodium falciparum, a causative agent of malaria, is a well-characterized obligate intracellular parasite known for its ability to remodel host cells, particularly erythrocytes, to successfully persist in the host environment. However, the current levels of understanding from the laboratory experiments on the host-parasite interactions and the strategies pursued by the parasite to remodel host erythrocytes are modest. Several computational means developed in the recent past to predict host-parasite/pathogen interactions have generated testable hypotheses on feasible protein-protein interactions. We demonstrate the utility of protein structure-based protocol in the recognition of potential interacting proteins across P. falciparum and host erythrocytes. In concert with the information on the expression and subcellular localization of host and parasite proteins, we have identified 208 biologically feasible interactions potentially brought about by 59 P. falciparum and 30 host erythrocyte proteins. For selected cases, we have evaluated the physicochemical viability of the predicted interactions in terms of surface complementarity, electrostatic complementarity, and interaction energies at protein interface regions. Such careful inspection of molecular and mechanistic details generates high confidence on the predicted host-parasite protein-protein interactions. The predicted host-parasite interactions generate many experimentally testable hypotheses that can contribute to the understanding of possible mechanisms undertaken by the parasite in host erythrocyte remodeling. Thus, the key protein players recognized in P. falciparum can be explored for their usefulness as targets for chemotherapeutic intervention.

Enriching the annotation of Mycobacterium tuberculosis H37Rv proteome using remote homology detection approaches: insights into structure and function.

The availability of the genome sequence of Mycobacterium tuberculosis H37Rv has encouraged determination of large numbers of protein structures and detailed definition of the biological information encoded therein; yet, the functions of many proteins in M. tuberculosis remain unknown. The emergence of multidrug resistant strains makes it a priority to exploit recent advances in homology recognition and structure prediction to re-analyse its gene products. Here we report the structural and functional characterization of gene products encoded in the M. tuberculosis genome, with the help of sensitive profile-based remote homology search and fold recognition algorithms resulting in an enhanced annotation of the proteome where 95% of the M. tuberculosis proteins were identified wholly or partly with information on structure or function. New information includes association of 244 proteins with 205 domain families and a separate set of new association of folds to 64 proteins. Extending structural information across uncharacterized protein families represented in the M. tuberculosis proteome, by determining superfamily relationships between families of known and unknown structures, has contributed to an enhancement in the knowledge of structural content. In retrospect, such superfamily relationships have facilitated recognition of probable structure and/or function for several uncharacterized protein families, eventually aiding recognition of probable functions for homologous proteins corresponding to such families. Gene products unique to mycobacteria for which no functions could be identified are 183. Of these 18 were determined to be M. tuberculosis specific. Such pathogen-specific proteins are speculated to harbour virulence factors required for pathogenesis. A re-annotated proteome of M. tuberculosis, with greater completeness of annotated proteins and domain assigned regions, provides a valuable basis for experimental endeavours designed to obtain a better understanding of pathogenesis and to accelerate the process of drug target discovery.

From workstations to workbenches: Towards predicting physicochemically viable protein-protein interactions across a host and a pathogen.

The understanding of protein-protein interactions is indispensable in comprehending most of the biological processes in a cell. Small-scale experiments as well as large-scale high-throughput techniques over the past few decades have facilitated identification and analysis of protein-protein interactions which form the basis of much of our knowledge on functional and regulatory aspects of proteins. However, such rich catalog of interaction data should be used with caution when establishing protein-protein interactions in silico, as the high-throughput datasets are prone to false positives. Numerous computational means developed to pursue genome-wide studies on protein-protein interactions at times overlook the mechanistic and molecular details, thus questioning the reliability of predicted protein-protein interactions. We review the development, advantages, and shortcomings of varied approaches and demonstrate that by providing a structural viewpoint in terms of shape complementarity and interaction energies at protein-protein interfaces coupled with information on expression and localization of proteins homologous to an interacting pair, it is possible to assess the credibility of predicted interactions in biological context. With a focus on human pathogen Mycobacterium tuberculosis H37Rv, we show that such scrupulous use of details at the molecular level can predict physicochemically viable protein-protein interactions across host and pathogen. Such predicted interactions have the potential to provide molecular basis of probable mechanisms of pathogenesis and hence open up ways to explore their usefulness as targets in the light of drug discovery.